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The measurement of the polarized gluon distribution function Ag(x) from pho- 
ton gluon fusion processes in electron proton deep inelastic scattering producing 
two jets has been investigated. The study is based on the MEPJET simulation 
program. The size of the expected spin asymmetry and corresponding statistical 
uncertainties for a possible measurement with polarized beams of electrons and 
protons at HERA have been estimated. The results show that the asymmetry can 
reach a few percents. 



1 Introduction 

After confirmation of the surprising EMC result, that quarks carry only a little 
fraction of the nucleon spin, this subject is actively Jpeing studied by several 
fixed target experiments at CERN, DESY and SLAC El Sofar only the polarized 
structure functions g\(x g ,Q 2 ) and 32(^9, Q 2 ) have been measured. These do 
not allow to distinguish between the role of quarks and gluons distributions. 
The measurement of the polarized gluon distribution Ag(x g ,Q 2 ) has become 
the key experiment in order to understand the QCD properties of the spin of 
the nucleon. We study here the possible direct measurement of Ag(x g ,Q 2 ) 
from jet events at the HERA collider, in the scenario where both the electron 
and the proton beam are polarized. 

The gluon distribution enters at leading order (LO) in the two-jets pro- 
duction cross section |] in deep inelastic scattering (DIS) (see Fig.[l]) , and the 
unpolarized gluon distribution g(x g ,Q 2 ) has indeed already been extracted 
from two-jets events by the HI collaboration at HERA. With the modest lu- 
minosity of 0.24 plr 1 obtained in 1993, first data on x g g(x g ) were extracted 
from dijets eventsB at LO, in a wide x g range 0.002 < x g < 0.2, at a mean Q 2 
of 30 GeV 2 . These results happen to be in good agreement with the gluon 
distribution extracted at LO from scaling violations of the structure function 
F 2 - 



a In the following the jet due to the beam remnant is not included in the number of jets. 
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Figure 1: Foynman diagrams for the dijet cross section at LO. 

2 Jet cross sections in DIS 

Deep inelastic electron proton scattering with several partons in the final state, 

e - (Z) +p(P) — > e~(l') + remnant(p r ) + parton l(pi) + . . . + parton n(p n ) (1) 

proceeds via the exchange of an intermediate vector boson V = 7, Z. In the 
following, Z-exchange will be neglected. We denote the momentum of the 
virtual photon, 7*, by q = 1 — 1', its absolute square by Q 2 , and use the 
standard scaling variables XBj = Q 2 / (2P ■ q) and y = P ■ q/P ■ I. The general 
structure of the unpolarized n-jet cross section in DIS is given by 

da had [n - jet] = ^ / dx a f a {x a ,H F ) da a {p = x a P, a 3 (fi 2 R ), H% fi 2 F ) (2) 

where the sum runs over incident partons a = q,q, g which carry a fraction 
x a of the proton momentum. a a denotes the partonic cross section from 
which collinear initial state singularities have been factorized out (in next-to- 
leading order (NLO)) at a scale /j,f and implicitly included in the scale depen- 
dent parton densities f a (x a , p? F ). For polarized lepton hadron scattering, the 
hadronic (n-jet) cross section is obtained from Eq. (|^) with the replacements 
(a had if a ,<r a ) — > (Aa had i A/ a , Aer°). The polarized hadronic cross section is 



2 



defined by Aa had = a^? — a^ d , where the left arrow in the subscript denotes 
the polarization of the incoming lepton with respect to the direction of its mo- 
mentum. The right arrow stands for the polarization of the proton parallel or 
anti-parallel to the polarization of the incoming lepton. The polarized parton 
distributions are defined by Af a (x a ,[i 2 F ) = f a ^(x a ,^ 2 F ) - f a {{x a ,^ F )- Here, 
fai(fai) denotes a parton a in the proton whose spin is aligned (anti-aligned) 
to the proton's spin. Aa a is the corresponding polarized partonic cross section. 

In Born approximation, the subprocesses 7* + q — > q + g, 7* + q — » q + g, 
7* +g — » q + q contribute to the two-jet cross section (Fig.^). The boson gluon 
fusion subprocess 7* + g — > q + q dominates the two jet cross section at low 
xbj (see below) and allows for a direct measurement of the gluon density in 
the proton. The full NLO corrections for Iwo jet production in unpolarized 
lepton hadron scattering are now available a and implemented in the ep — > n 
jets event generator MEPJET, which allows to analyze arbitrary jet definition 
schemes and general cuts in terms of parton 4-momenta. 

First discussions about jet production in polarized lepton-hadron scatter- 
ing can be found in Ref. [4], where the jets were defined in a modified "JADE" 
scheme. However, it was foundta that the theoretical uncertainties of the two- 
jet cross section for the "JADE" scheme can be very large due to higher order 
effects. These uncertainties are small for the cone scheme and the following 
results are therefore based on the cone algorithm, which is defined in the labo- 
ratory frame. In this algorithm the distance AR — J (A77) 2 + (A<fr) 2 between 
two partons decides whether they should be recombined to a single jet. Here 
the variables are the pseudo-rapidity 77 and the azimuthal angle (f>. We recom- 
bine partons with AR < 1. Furthermore, a cut on the jet transverse momenta 
°f Pr(j) > 5 GeV in the lab frame and in the Breit frame is imposed. We 
employ the one loop (two loop) formula for the strong coupling constant in a 
LO (NLO) analyses with a value for A^gL consistent with the value from the 
parton distribution functions. In addition a minimal set of general kinematical 
cuts is imposed on the virtual photon and on the final state electron and jets. 
If not stated otherwise, we require 5 GeV 2 < Q 2 < 2500 GeV 2 , 0.3 < y < 1, 
an energy cut of E(e') > 5 GeV on the scattered electron, and a cut on the 
pseudo-rapidity 77 — — lntan(0/2) of the scattered lepton (jets) of \t]\ < 3.5 
(|ry| < 2.8). These cuts are compatible with the existing detectors HI and 
ZEUS. 

Let us briefly discuss the choice of the renormalization and factorization 
scales nn and \ip in Eq. (Q). Jet production is a multi-scale problem and it is 
not a priori clear at,which scale a s and the parton densities are probed. How- 
ever, it was arguedBB, that the "natural" scale for jet production in DIS is set 
by the average k® of the jets in the Breit frame. Here, {k F (j)) 2 is defined by 
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2E?(1 — cosOjp), where the subscripts j and P denote the jet and proton, re- 
spectively (all quantities are defined in the Breit frame). It can be shownlj that 
J2j ^t(j) smoothly interpolates between the correct limiting scale choices, it 
approaches Q in the parton model limit and it corresponds to the jet trans- 
verse momentum (with respect to the 7*-proton direction) when the photon 
virtuality becomes negligible. It therefore appears to be the "natural" scale 
for multi jet production in DIS and we choose n 2 R = /if, = 1/4 (J^ . k^(j)) 2 in 
Eq. (||) as our standard choice. 

Let us first discuss some results for unpolarized dijet cross sections. If not 
stated otherwise, the lepton and hadron beam energies are 27.5 and 820XJeV, 
respectively. With the previous parameters and GRV parton densities one 
obtains a LO (NLO) two jet cross section a had (2 - jet) of 1515 pb (1470 pb). 
Thus the higher order corrections are small. This is essentially due to the rel- 
atively large cuts on the transverse momenta of the jets. As mentioned before, 
the boson-gluon fusion subprocess dominates the cross section and contributes 
to 80% to the LO cross section. 

In order to investigate the feasibility of the parton density determination, 
Fig. 2a shows the Bjorken Xgj distribution of the unpolarized two jet exclusive 
cross section. The gluon initiated subprocess clearly dominates the Compton 
process for small xgj in the LO predictions. The effective if -factor close to 
unity for the total exclusive dijet cross section is a consequence of compensating 
effects in the low x (K > 1) and high x (K <1) regime. 

For the isolation of parton structure functions we are interested in the 
fractional momentum x a of incoming parton a (a = g, g), however, and in 
dijet production xbj and x a differ substantially. Denoting as Sij the invariant 
mass squared of the produced dijet system, and considering two-jet exclusive 

events only, the two are related by x a — XBj ^1 + ■ The s.y distribution 

of Fig. 2b exhibits rather large NLO corrections as well. The invariant mass 
squared of the two jets is larger at NLO than at LO (the mean value of se- 
nsing to 620 GeV 2 at NLO from 500 GeV 2 at LO). 

The NLO corrections to the XBj and Sy distributions have a compensating 
effect on the x a distribution in Fig. 2c, which shows similar shapes at LO and 
NLO. At LO a direct determination of the gluon density is possible from this 
distribution, after subtraction of the calculated Compton subprocess. This 
simple picture is modified in NLO, however, and the effects of Altarclli-Parisi 
splitting and low px partons need to be taken into account more carefully to 
determine the structure functions at a well defined factorization scale /if in 
NLO. 

In the following we discuss some results for polarized dijet production. Our 
standard set of polarized parton distributions is "gluon, set A" of Gehrmann 
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and StirlingEI, for which J Ag = 1.8 at Q 2 = 4 GeV 2 . Using the same kinemat- 
ical cuts as before, the LO polarized dijet cross section Acr(2 — jet) is -45 pb. 
This negative result for the polarized dijet cross section is entirely due to the 
boson-gluon fusion process, which is negative for xgj < 0.025. The contribu- 
tion to the total polarized dijet cross section is -53 pb. The contribution from 
the quark initiated subprocess is positive over the whole kinematical range and 
contributes with 8 pb to the resulting dijet cross section. Note, however, that 
the shape of the x g distribution in the polarized gluon density is hardly (or 
even not at all) constrained by currently available DIS data (in particular for 
small x g ). Alternative parametrizations of the polarized gluon distributions 
in the small x g region, which are still consistent with all present datacl, can 
lead to polarized cross-sections Aa had , which are at least a factor two larger 
than those obtained with "gluon, set A" . It is clear from Fig. 2d, that the dijet 
events are particularly sensitive to this lower x g range. Note, that this frac- 
tional momentum distribution is again related to xgj by x a — xgj I 1 + J . 



The corresponding xsj and Sy distributions are not shown here. 
3 Experimental asymmetries 

In order to study the feasibility and the sensitivity of the measurement of the 
spin asymmetry at HERA, we have assumed polarizations of 70% for both the 
electron and the proton beams and statistical errors were calculated for a total 
luminosity of 200 pb^ 1 . 

The expected experimental asymmetry < A >— A lad^'et^ unc ^ er these 
conditions is shown in Fig. 3 EL cLS cl function of xsj and in Fig. ||b-d as a 
function of x g . Figures a) and b) correspond to the nominal kinematical cuts 
defined previously, except for the Q 2 range which was extended to lower values 
2 < Q 2 < 2500 GeV 2 . The cross section integrated over all variables is 2140 pb, 
where 82% of the contribution comes from gluon-initiated events. The asym- 
metry averaged over all variables is < A >— —0.015 ± 0.0015. It is negative 
at low XBj and becomes positive at XBj > 0.01. 

In figures c) and d), a further cut was made: Q 2 < 100 GeV 2 . This cut 
permits to reject the positive contributions to the asymmetry coming from high 
Q 2 (equivalent to high xsj ) events, where the contribution of quark-initiated 
events is higher. All the remaining events were separated in two bins in s,y 
-the invariant mass of the dijet- and two bins in y, as the asymmetry is very 
sensitive to these two variables. Fig. c) corresponds to low invariant masses 
(sij < 500 GeV 2 ), and Fig. d) to high ones (sy > 500 GeV 2 ). Open points 
show low y values ( y < 0.6), and closed points, high y values ( y > 0.6). In 
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l°S.o(*a> log,„(x a ) 



Figure 2: a) Dependence of the unpolarized two-jet cross section on Bjorken XBj for the 
quark and gluon initiated subprocesses and for the sum. Both LO (dashed) and NLO (solid) 
results arc shown; b) Dijet invariant mass distribution in LO (dashed) and in NLO (solid) 
for unpolarized dijet production; c) Same as a) for the x a distribution, x a representing 
the momentum fraction of the incident parton at LO; d) Dependence of the LO polarized 
two-jet cross section on x a for the quark and gluon initiated subprocesses (dashed) and for 

the sum (solid). 

the best case, the asymmetry reaches values as high as 12% (Fig. d). 

Reducing the beam energy to 410 GeV, instead of the nominal 820 GeV, 
does not improve the signal in average, although the mean value of y is higher. 
The asymmetry signal increases only for a few points around x g > 0.1, since a 
higher incident energy probes slightly higher values of x g . 

The results show that if the assumed luminosity and beam polarizations 
can be delivered at HERA, the present detectors HI and ZEUS will be in a 
comfortable position to measure a spin asymmetry of a few per cent in average, 
with a few per mil statistical precision. In order to minimize the experimen- 
tal systematic uncertainties, it is desirable to have in the HERA ring bunch 
trains of protons with alternate helicity. On the theoretical side, NLO QCD 
corrections are needed. The NLO corrections reduce the renormalization 
and factorization scale \i f dependence (due to the initial state collinear factor- 
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Figure 3: Expected asymmetries as a function of xsj (a) and x g (b-d) for a luminosity of 
200 pb^ 1 and beam polarizations P e = P P = 70%; Fig. a) and b):nominal cuts; Fig.c) and 
d): Q 2 < 100 GeV 2 ; In Fig.c) and d) data are separated in bins of Sij and y. 
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ization) in the LO calculations and thus reliable predictions in terms of a well 
defined strong coupling constant and scale dependent parton distributions be- 
come possible. At the moment, these corrections are only available for unpolar- 



ized jet production BH. One expects for the asymmetry < A > 



_ Ao-'"" i (2-jet) 
T had (2-jet) ' 

that the scale dependence in the individual cross sections partly cancels in the 
ratio. In fact, varying the renormalization and factorization scales between 
li\ = [i} F = 1/16 (Ej ^t(j)) 2 and n\ = m| = 4 (E, H)? m the LO cross 
sections introduces an uncertainty for the ratio A of less than 2 %, whereas 
the uncertainty in the individual cross sections is much larger. 

In conclusion, the dijets events from polarized electron and polarized pro- 
ton collisions at HERA can provide a good measurement of the gluon polariza- 
tion distribution for 0.002 < x g < 0.2, the region where x g Ag(x g ) is expected 
to show a maximum. 
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